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ABSTRACT 

The symbiotic X-ray binary (SyXB) 4U 1954+319 is a rare system hosting a peculiar neutron star (NS) and an 
M-type optical companion. Its ~5.4 hr NS spin period is the longest among all known accretion-powered pulsars 
and exhibited large (~7%) fluctuations over 8 yr. A spin trend transition was detected with Swift /BAT around 
an X-ray brightening in 2012. The source was in quiescent and bright states before and after this outburst based 
on 60 ks Suzaku observations in 2011 and 2012. The observed continuum is well described by a Comptonized 
model with the addition of a narrow 6.4 keV Fe-Ka line during the outburst. Spectral similarities to slowly rotating 
pulsars in high-mass X-ray binaries, its high pulsed fraction (~60%-80%), and the location in the Corbet diagram 
favor high 5-field (>10 12 G) over a weak field as in low-mass X-ray binaries. The observed low X-ray luminosity 
(10 33 -10 35 ergs -1 ), probable wide orbit, and a slow stellar wind of this SyXB make quasi- spherical accretion in 
the subsonic settling regime a plausible model. Assuming a ~10 13 G NS, this scheme can explain the ~5.4 hr 
equilibrium rotation without employing the magnetar-like field (~10 16 G) required in the disk accretion case. The 
timescales of multiple irregular flares (~50 s) can also be attributed to the free-fall time from the Alfven shell for 
a ~10 13 G field. A physical interpretation of SyXBs beyond the canonical binary classifications is discussed. 
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1. INTRODUCTION 

Symbiotic X-ray binaries (SyXBs) are a growing new sub- 
class of X-ray binaries which consist of an X-ray bright neutron 
star (NS) and an M-type giant primary star. Originally SyXBs 
were seen as a subclass of low-mass X-ray binaries (LMXBs), 
where accretion typically is from a K-type optical counterpart 
with a mass M c < 1 M©. It took more than 30 yr, however, to 
recognize that their properties are different from LMXBs and 
also from those of high-mass binaries (HXMBs), i.e., systems 
where the mass donor is an O- or B-type giant with M c > 10 M 0 . 

SyXBs are now considered as a separate subclass of X-ray bi- 
naries. They are named in analogy to symbiotic binaries, which 
consist of a white dwarf and an M-giant companion, but their 
nature has not yet been understood. 

Some SyXBs were known as famous X-ray sources since 
the early X-ray observations: GX 1+4 (Davidsen et al. 

1977; Makishima et al. 1988; Chakrabarty & Roche 1997), 

4U 1700+24 (Garcia et al. 1983; Masetti et al. 2002), and Set 
X-l (Kaplan et al. 2007). With the discovery of new sources 
such as IGR J16 194— 2810, there are now nearly five SyXBs 
and some candidates in our Galaxy; e.g., five objects out of 
ten candidates (Table 1 in Lii et al. 2012) have been confirmed 
as SyXBs, while three were ruled out (see, e.g., Masetti et al. 

2012). Observational features of the SyXB class are (1) quite 
long orbital periods, 8 * e.g., 1161 days in GX 1+4; (2) long NS 


8 For many SyXB, the orbital period has not yet been known nor confirmed 

in long-term X-ray studies (e.g., Corbet et al. 2008). 


spin periods over ~ 110-18300 s as shown in Figure 1; and 
(3) high X-ray variability ranging from short to long timescales 
(~1 s to up to a year). 

SyXB X-ray spectra are in general well described by an 
absorbed power law with a photon index T ~ 0.5-2. 0 and 
a high energy roll-over, which is often modified in the soft 
X-rays by strong absorption. This spectrum resembles that 
of X-ray pulsars in HMXBs, and thus, SyXBs have been 
interpreted in the literature as NSs with ~10 12 G field in 
HMXBs. As another interpretation, following the classical 
LMXB classification based on optical companions, SyXBs were 
also interpreted as NSs in the low-luminosity low/hard state of 
LMXB as having weak magnetic field (5-field), e.g., <10 10 G 
(Nagae et al. 2008; Kitamura et al. 2014). In any of these cases, 
the NSs in SyXB are expected to have 5-fields strengths of 
around 10 12 G or much lower (e.g., <10 10 G). 

The canonical 5-fields implied by the spectral analogy with 
HMXBs and LMXBs, however, are in contradiction to the fields 
implied by their timing behavior. The slow NS pulsations in 
SyXBs are sometimes in spin equilibrium or show a large spin- 
up rate (Makishima et al. 1988; Gonzalez-Galan et al. 2012). 
If interpreted in terms of standard disk accretion torque theory 
(Ghosh & Lamb 1979), the slow rotation and torque reversals 
imply extraordinary strong dipole fields (~ 10 14-16 G), which are 
close to or exceed the quantum critical field, 5 cr = m 2 c 3 /(eh) = 
4.4 x 10 13 G. Such a magnetar scenario for binary systems has 
been proposed, e.g., to explain the 1.6 hr pulsation and the 
high spin-up rate of the SyXB IGR J16358 — 4726 (Patel et al. 
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Figure 1 . Distribution of NS rotation periods in HMXBs (Liu et al. 2006), 
LMXBs (Liu et al. 2007), and confirmed SyXBs. Measured rotation periods of 
SyXBs are -18,300 s (4U 1954+319), -110 s (Set X-l), and -150 s (GX 1+4). 
The 4U 1954+319 is in the HMXB range. 

(A color version of this figure is available in the online journal.) 


2007). Similar scenarios have also been discussed for slowly 
rotating X-ray pulsars in HMXBs such as 4U 0114+65 (Li & 
van den Heuvel 1999) or 4U 2206+54 (Reig et al. 2012), highly 
variable supergiant fast X-ray transients (Bozzo et al. 2008), and 
some slowly rotating NSs in the Small Magellanic Cloud (Ho 
et al. 2014; Klus et al. 2014). Except for the timing properties, 
however, no direct evidence for ultra-strong 5-fields has been 
found so far (e.g., Ikhsanov 2007). 

As shown by Shakura et al. (2012) and Postnov et al. (2010, 
2012), the extremely strong fields can be avoided when one as- 
sumes that the accretion does not happen through an accretion 
disk but via quasi- spherical accretion. In these models, matter 
is gravitationally captured, e.g., from the donor’s stellar wind. 
For low mass accretion rates, M < 4 x 10 16 g s _1 , gravitation- 
ally captured matter subsonically accretes onto the NS forming 
an extended, quasi-spherical shell around the magnetosphere, 
where large scale convective motions and subsonic turbulence 
lead to accretion onto the NS. Contrary to the canonical high- 
luminosity Bondi-Hoyle accretion, where matter supersonically 
accretes and the sign of angular momentum (prograde or ret- 
rograde) determines the spin behavior (Perna et al. 2006), for 
quasi-spherical accretion the spin-up/ spin-down is determined 
by the specific angular momentum of matter at the magneto- 
spheric boundary and the angular velocity of the NS. 

In this paper, we discuss Suzaku observations of the SyXB 
4U 1954+319 hosting the slowest rotating NS in an X-ray 
binary system. This source was originally discovered by Uhuru 
(Forman et al. 1978) and observed in the 1980s with Ariel 
(Warwick et al. 1981), EXOSAT , and Ginga (Tweedy et al. 
1989). These observations suggested that the source was a wind 
accreting HMXB with a very inhomogeneous wind (Tweedy 
et al. 1989), although no optical counterpart was known at that 
time. The source had virtually been forgotten until, after 20 yr of 
no observations, Chandra' s localization of the source position 
allowed the optical identification of the donor, which was found 
to be an M4-5 III star (M c ~ 1.2 M 0 ) at a distance of 1.7 kpc 


(Masetti et al. 2006b). Therefore, 4U 1954+319 was re-classified 
as a SyXB. 

Further progress in our understanding of the characteristics 
came with the discovery of a ~5 hr periodicity in Swift /BAT 
monitoring data by Corbet et al. (2006, 2008). This period 
is strongly variable (Marcu et al. 2011). So far, no orbital 
period has been discovered, although Mattana et al. (2006) argue 
that the lower limit of the orbital period is ~400 days. These 
parameters are incompatible with typical white dwarf systems, 
such that the compact object in 4U 1954+319 is likely a NS. 
As shown in Figure 1, the extremely long spin period makes 
4U 1 954+3 1 9 the slowest rotating accretion-powered NS system 
known to date 9 and a prominent example for understanding the 
SyXB class. 

Currently, no consensus has yet been obtained on the accretion 
models for SyXBs, or the 5-field of their NSs. Revisiting 
and developing the quasi-spherical model further requires the 
spectral and timing investigations of SyXBs. In this paper, we 
report on the broadband timing and spectral properties of 4U 
1954+319 and discuss a possible accretion model beyond the 
conventional HMXB and LMXB classification. 

2. OBSERVATION AND DATA REDUCTION 
2.7. Suzaku Observations 

There are two pointed Suzaku (Mitsuda et al. 2007) observa- 
tions of 4U 1954+319. The first one, in A06, was conducted 
on 2011 October 23-24 for a gross on-source time of 117 ks. 
Nearly 1 yr later, on 2012 October 5, MAXI monitoring detected 
a brightening of this source, 10 with the X-ray intensity reaching 
~15 mCrab in the 3-10 keV band. About 27 days later, on 2012 
November 1/2, a follow-up target of opportunity (ToO) obser- 
vation was performed with Suzaku for a gross on-source time of 
120 ks. Table 1 summarizes the two Suzaku observations. 

The X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
on board Suzaku , coupled to the X-ray Telescopes (XRT; 
Serlemitsos et al. 2007), consists of two front-illuminated (FI) 
CCD detectors, XIS0 and XIS3, and one back-illuminated 
sensor (XIS1). The CCDs are sensitive in the 0.2-12 keV range 
and, when operated in the full window mode, read out every 
8 s. We also used data from the hard X-ray detector (HXD; 
Takahashi et al. 2007), which consists of Silicon PIN diodes 
(HXD-PIN) and GSO scintillators (HXD-GSO) and covers the 
10-70 keV and 50-600 keV bands, respectively. The HXD was 
operated in the nominal mode during both observations with 
61 /xs time resolution. The bright X-ray source Cygnus X- 
1, which is only 3? 15 from 4U 1954+319, is still outside the 
HXD-PIN ’s tightly collimated field of view (FOV; 34' x 34' 
FWHM) but is within the FOV of the HXD-GSO (4° x 4°). 
Thus, we did not use the HXD-GSO data in the following 
analyses due to a possible contamination from Cyg X-l. Based 
on the INTEGRAL source catalog, there is no other strong 
contamination source within the HXD-PIN FOV. 

2.2. Data Processing 

We analyzed the XIS and HXD data sets as obtained from 
the standard pipeline processing, revision 2.7 and 2.8 for the 


9 A long periodicity, P ~ 6.67 hr, was also detected from the central object 
1E161348— 5055 in the young shell-type supernova remnant RCW 103. This 
object, however, has not been confirmed as an accretion-powered pulsar. 

10 Trigger ID 6206800002: 

http://maxi.riken.jp/pipermail/x-ray-star/2012-0ctober/000187.html. 
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Table 1 

Log of Suzaku Observations of 4U 1954+319 



2011 

A06 

2012 
A07 ToO 

Observation start 
Start MJD 
Observation end 
End MJD 
Nominal position 
ObsID 

Oct. 23 08:05 
55857.337 
Oct. 24 17:17 
55858.720 
XIS 

406046010 

Nov. 1 02:39 
56232.110 
Nov. 2 12:02 
56233.501 
XIS 

907005010 

XIS signal rate (cnt s -1 ) 
XIS exposure (ks) 

PIN signal rate (cnt s -1 ) 
PIN exposure (ks) 

X-ray Intensity 

0.18/0.18 

60.2 

53.7 

0.28 mCrab 

5. 1/4.5 
60.8 
0.70 
50.5 

9. 5 mCrab 

F x (10 -12 erg s -1 cm -2 ) 2-10 keV 
F x (10“ 12 erg s“‘ cm” 2 ) 1-70 keV 

5.36 ±0.03 

6-28: 2 o 6 5 9 2 

190.0^° 0 2 

600.4+ 3 4 '\ 

L x (erg s” 1 ) 1-70 keV 

2.2 x 10 33 

2.1 x 10 35 


Notes. The 0.5-10 keV XIS and 12-70 keV HXD-PIN signal rates are averaged 
during the observation. The XIS rates are shown for the XIS-FI and -BI CCD 
imagers. The source and background regions are stated in Section 2.2. The 
X-ray intensity is evaluated in the 2-10 keV band. F x and L x are evaluated 
using the comptt model. 

2011 and 2012 data sets, respectively. All data analyses were 
performed using HE ADAS version 6.1 1 or later. Following the 
standard Suzaku screening criteria, X-ray events were discarded 
if one of the following conditions held true: (1) they were 
recorded during the South Atlantic Anomaly, TSAA < 436 s 
and TSAA ^ 500 s for the XIS and the HXD, respectively, and 
later than 180 s before the next passage for the HXD; (2) while 
the satellite was in regions of low geomagnetic cutoff rigidity 
(^6 GV); (3) while the target was closer than 5° to Earth’s rim; 
(4) if the instantaneous pointing direction deviated more than 
1 '.5 from the mean; or (5) during periods of telemetry saturation. 
After screening, the net exposures for the 2011 observation were 
60.2 ks (XIS) and 53.7 ks (HXD), while those in 2012 were 
60.8 ks (XIS) and 50.5 ks (HXD), respectively. 

On- source XIS events were extracted from a region of 
2^5 radius centered on the source. Background events were 
derived from an annulus 4^0-7^5 around the source. From the 
HXD-PIN data, we subtracted the non X-ray Background 
(NXB) created with the LCFITD method (Fukazawa et al. 2009) 
and filtered following the same screening criteria as those used 
for the on-source data. The cosmic X-ray Background (CXB) 
was modeled as described by Enoto et al. (2010a) with the same 
spectral shape as given by Moretti et al. (2009). In the following 
analyses, unless otherwise specified, all uncertainties are given 
at the 68% (la) confidence level. 

2.3. Reduction of the X-Ray Spectra 

For the spectral analysis, we summed the two XIS-FI sensors 
and utilized XIS data from the 0.8-10 keV band and the 
HXD-PIN data from the 15-70 keV band. The 1 .7-1 .9 keV XIS 
data are discarded due to calibration uncertainties. Uncertainties 
in the flux normalization of the instruments are taken into 
account with a multiplicative constant. The cross-normalization 
factor of HXD-PIN relative to XIS-FI was fixed at 1.171 
based on the XlS-nominal Crab Nebula calibration, 11 while 


1 1 http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/ 
suzakumemo-2008-06.pdf 


that of XIS-FI and XIS-BI was allowed to differ by up to 
1%. A 0.5% systematic uncertainty was added in quadrature 
to the 2012 XIS spectral bins in order to absorb instrumental 
calibration uncertainties. The response matrix files and auxiliary 
response file for the XIS were produced using the FTOOFs 
xisrmfgen and xissimarfgen (Ishisaki et al. 2007), while 
standard detector response files (epoch 11) were used for the 
HXD-PIN. 

3. SUZAKU OBSERVATIONS OF 4U 1954+319 
3.1. Long-term Variability 

Figure 2 shows the long-term variability of the X-ray intensity 
of 4U 1954+319, covering ~17 yr since 1994. X-ray count rates 
were re-binned to 6 and 12 days resolution and converted into 
the mCrab unit using conversion factors provided by the RXTE/ 
All-Sky Monitor (ASM) and Swift/ Burst Alert Telescope (BAT) 
teams. 12 As shown in the 20 day running average of the fluxes 
(Figure 2, blue curve), the X-ray intensity of 4U 1954+319 is 
highly variable on a timescale of years, with fluxes changing 
from a few mCrab to at least 20 mCrab or more. In addition, 
faster variability on timescales of weeks or days is also seen 
with irregular flares. 

Figure 2 also shows the spin history of 4U 1954+319 as 
measured with Swift /BAT. We determined the pulse periods 
using standard epoch folding (Feahy et al. 1983) and using 
pulse profiles with 19 phase bins. Spin values are shown with 
filled circles if pulses are detected with x^d > 3- The measured 
period history is consistent with previous work covering data 
taken before 2010 (Marcu et al. 2011). The present work extends 
the analysis until the end of 2012 and shows the detection 
of a new change from spin-down to spin-up which occurred 
around 2012 March (MJD ~ 56000), when the average X-ray 
intensity increased to ~25 mCrab. The analysis shows that the 
first and second Suzaku observations were taken during the spin- 
down and spin-up phases, respectively. The former observation 
was almost in quiescence (~0.3 mCrab), before the outburst, 
while the latter was during a brighter phase after the flaring 
(~10 mCrab) and the spin-transition. 

3.2. Timing Analysis 

Figure 3 shows the XIS and HXD-PIN light curves of 4U 
1954+319 from the two observations. Similar to the behavior 
observed on long timescales (Figure 2), the X-ray light curve 
is highly variable also on short timescales. Figure 4 further 
illustrates the enlarged light curve. The X-ray fluctuation was 
moderate in 2011, while spiky flare events were seen in 2012 
with <100 s scale, showing a rapid flux change by an order of 
magnitude or more. During the 2012 outburst (Figure 4), the 
spiky flares become enhanced in the latter part of observation. 
An increase of the short-term fluctuation is also seen in the 
power spectral density (PSD) in Figure 5. The PSD continuum 
can be described by a power-law components with best-fit slopes 
of —0.77 in 2011 and —1.66 in 2012. There are no signatures 
of clear short-term periods or quasi periodic oscillations. No 
low-frequency roll over was detected, either. 


12 The long-term RXTE / ASM monitoring data are provided by the 
RXTE / ASM team (http://xte.mit.edu), while the Swift /BAT transient 
monitor results provided by the Swift/BAT team 
(http://heasarc.gsfc.nasa.gov/docs/swift/results/transients/). The reported 
conversion factors, assuming a Crab-like spectrum, are 
7.5 x 10 -2 counts s -1 mCrab -1 and 2.2 x 10 -4 counts cm -2 s -1 mCrab -1 , for 
ASM and BAT, respectively. 
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Figure 2. Long-term light curve of 4U 1954+319. (a) 2-10 keV RXTE / ASM (triangle, green) and 15-50 keV Swift /BAT (filled circle, orange) light curves of 4U 
1954+319 since 1994. The solid blue curve represents a running average of the X-ray intensity with a window size of 20 days. The detected pulse period history by 
the Swift /BAT is overlaid as black filled circles (right hand axis). Two MAXI detections of the flaring activity of this source are also shown: MAXI trig 1 on 2012 
June 29 (Sugizaki et al. 2010) and trig 2 on 2012 October 5 (see the text). The two Suzaku observations are indicated as red star symbols, (b) Enlarged version of the 
light curves around the two Suzaku observations, showing the 7 days averaged X-ray intensity from Swift /BAT (orange). 

(A color version of this figure is available in the online journal.) 



Figure 3. Background- subtracted light curves of 4U 1954+319 obtained with Suzaku in 2011 October (left) and 2012 November (right). Gray strips represent the non 
good time intervals (non-GTIs) during the observation. From top to bottom, the panels show the 0.5-1. 0 keV XIS (binning 360 s (2011) and 120 s (2012)), 1-10 keV 
XIS (64 s and 32 s), and 15-70 keV HXD-PIN (160 s and 80 s) count rates. The XIS background was subtracted as described in Section 2.2, while the simulated NXB 
(~0.28 count s -1 ) and CXB (~0.02 counts s -1 ) were subtracted from the dead-time-corrected HXD-PIN data. Estimated pulse shapes are overlaid in each panel, using 
the Swift /BAT period described in Section 3.2. 

(A color version of this figure is available in the online journal.) 


Epoch folding of the 2012 XIS and HXD data yielded a spin 
period of 5.68 ± 0.28 hr. However, the period determination 
with Suzaku is not accurate enough, due to the limited exposure 
which corresponds to only ~6 cycles of the NS rotation. Thus, 
we searched the Swift/ BAT monitoring data (Section 3.1) for 
the spin period using intervals of ±60 days and ±400 days 
around the Suzaku pointings. The pulse periods derived from 
these data are P = 5.70 ± 0.01 hr ( x 2 /v = 104.7/14 = 7.5) 
and P = 5.76 ± 0.01 hr (x 2 /v = 72.3/14 = 5.2), in 2011 and 
2012, respectively, where 15 phase bins were used for the epoch 
folding. Average pulse profiles, obtained at these periods, are 
superposed in Figure 3 on the Suzaku light curves. Even though 
sporadic flares are rather strong, we can trace individual pulses 
in most cases. 

Figure 6 shows the pulse profiles of the two Suzaku obser- 
vations. The pulse profile in 2012 is approximately sinusoidal 
below 1 keV, with a single peak at phase ~0.4. A second peak 
at phase ~0.2 appears at higher energies. Figure 7 shows the 
energy dependence of the pulsed fraction (PF) of 4U 1954+319, 


defined as PF — (^max ^min)/(^max + ^min)? where Rmax and 
R m i n are the maximum and minimum of the pulse profile. The 
fraction is 60%-80% in the 0.5-30 keV band and increases with 
energy. This is a typical trend often seen in X-ray pulsars in 
HMXBs (Futovinov & Tsygankov 2008). 

3.3. X-Ray Spectrum During the Flare in 2012 

Figure 8(a) shows Suzaku spectra of 4U 1954+319 obtained 
on the two occasions. We superpose the modeled HXD-PIN 
NXB and CXB, together with the typical 1% systematic uncer- 
tainty of the NXB. The HXD-PIN detection is marginal in 201 1. 
The ratio between the two spectra is shown in Figure 8(b). In 
2012, 4U 1954+319 was brighter by nearly 1-2 orders of mag- 
nitude and exhibited harder spectra. 

We now turn to modeling the X-ray spectrum in Figure 9, 
and summarize the best-fit parameters in Table 2. Fitting the 
data with a single, absorbed power law leaves a soft excess 
below 2 keV, an iron line at ~6.4 keV, and a spectral roll 
over (Figure 9(a)). Here we use the latest interstellar absorption 
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Figure 4. XIS 1-10 keV light curves during the 2011 (a, b) and the 2012 
observation (c, d). Data are binned to 8 s resolution. Time durations are 
(a) 22.0-25.8 ks, (b) 55.8-59.6 ks, (c) 13.6-17.4 ks, and (d) 88.2-92.0 ks 
in Figure 3. 

(A color version of this figure is available in the online journal.) 

model tbnew (Wilms et al. 2000) with an improved abundance 
wilm and cross section tables vern (Verner et al. 1996). 
Due to higher data quality compared to earlier INTEGRAL 
measurements, the broken power-law model used by Fiirst 
et al. (2011b) and Marcu et al. (2011) was not successful 
(Xy = 2373.6/697 = 3.40), either. The exponentially cutoff 
power law with an additional soft blackbody and an Fe-line, 
improved the fit to x 2 / v = 971.5/696 = 1.40 (Figure 9(b)) 
with a photon index of T = 1 .21 , a cutoff energy E cut = 20 keV, 
and a blackbody temperature of kT^ = 0.13 keV. The spectral 
parameters are similar to those found by Mattana et al. (2006) 
by applying the same model to the BeppoSAX data, although 
our fit is not satisfactory. 

As alternative representations of the roll-over of the contin- 
uum, we employed two empirical models: (1) a high energy 
cutoff model of the form (in XSPEC notation) 

phabs x (bbodyr ad + gauss + power law x highecut), (1) 

where 

1 for E < E c 

highecut - J exp[(£c _ E)/Ef] for E ^ Ec (2) 

with the cutoff energy E c ~ 4.2 keV and ^-folding energy 
Ef ~ 21.0 keV. We also modeled the data using the (2) a 
negative and positive power law with an exponential (NPEX), 

phabs x (bbodyrad + gauss + cutoff pi + cutof f pi), (3) 

which is often utilized to represent broadband continuum from 
the accretion column of highly magnetized NSs (Mihara 1995; 
Makishima et al. 1999; Enoto et al. 2008). The cutoff energy 
and the photon index of the NPEX continuum were obtained as 
6.6 keV and 0.77, respectively, and the blackbody temperature 
as 0.14 keV. As shown in Figures 9(c), (d), and Table 2, both 



Figure 5. Power spectrum densities (PSDs) of the 0.8-10 keV XIS data in 
2011 and 2012 observations as calculated with the FTOOL powspec using the 
background- subtracted 8 s binned light curve. The PSDs are normalized such 
that the white noise level corresponds to 2. The best-fit power-law models are 
shown with slopes of —0.77 and —1.66 in 2011 and 2012, respectively, when 
fitted above the white noise level. 

(A color version of this figure is available in the online journal.) 

models reproduced the broadband spectra quite well (x 2 ~ 1 .2). 
If we exclude the soft blackbody, x 2 becomes worse (x 2 > 3.5). 
Replacing the blackbody component with a diffuse plasma 
emission model (apec 13 ; Smith et al. 2001), did not give 
acceptable fits (x 2 > 1.9). 

Finally, as a more physically based interpretation, we 
also tried an optically thick Comptonization model, compTT 
(Titarchuk & Lyubarskij 1995). This gave an acceptable fit, 
X 2 = 780.8/697 = 1.12, with a null hypothesis probability of 
0.02 (Figure 9(e)). We therefore adopt this model as a canonical 
continuum model of 4U 1954+319. The derived best-fit param- 
eters are a soft photon (Wien) temperature of kTo ~ 1 keV, 
plasma temperature kT e ~ 8 keV, and an optical depth of 
r ~ 9.0. These derived kTo and kT e values are within the typical 
range of the comptt model of this source, kTo ~ 0.6-1. 3 keV 
and kT e ~ 3.0-13 keV, by the multi X-ray satellite studies 
(Masetti et al. 2007b) and by the INTEGRAL observation (Marcu 
et al. 201 1). Combining the previous two successful models, the 
broadband spectrum is explained by the Comptonized model. 
Figure 10 gives a vF v representation of the data and the model. 
The 2-10 keV and 1-70 keV absorbed X-ray fluxes are (1.90 zb 
0.01) x 10 -10 and (6.0 zb 0.1) x 1 0“ 10 erg s -1 cm -2 , translating 
to absorption-corrected X-ray fluxes of (2.03 =b 0.01) x 10“ 10 
and (6.3 =b 0.1) x 10“ 10 erg s -1 cm -2 , respectively. 

3.4. Iron Line and Search for CRSF 

The 6.4 keV neutral Fe-Ka line emission is obvious in 
Figure 9(a). In the compTT continuum fits, its center energy 
was found at 6.375 + qqq° 8 keV with an equivalent width (EW) 
of EW = 28.6 =b 0.3 eV and a total line flux of / Fe = 
7.5^0 7 x 10 -5 photon cm -2 s -1 . The line is narrow with a 
width a < 50 eV (la), corresponding to a Doppler width 
Voo ^ 2400kms _1 . 

Detailed inspection of the spectrum indicates a shallow dent 
at ~7. 1 keV, and if we let the iron abundance free, the data favor 
the 2.4 times larger iron abundance relative to the default im- 
proving the fit by Ax 2 = 15.3 for A v = —1 (F-test probability 


13 Astrophysical Plasma Emission Code, a collisonally ionzied diffuse gas 
model http://atomdb.org/. 
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Figure 6. Pulse profiles of 4U 1954+319 folded with 5.70 hr and 5.76 hr in 2011 (left) and 2012 (right), respectively. Panels (a), (b), (c), and (d) represent the 
0.5-1. 0 keV XIS, the 1-5 keV XIS, the 5-10 keV XIS, and the 15-70 keV HXD-PIN energy bands, respectively. Backgrounds were subtracted in the same way as for 
Figure 3. 

(A color version of this figure is available in the online journal.) 


Table 2 

Comparison of the 2012 and 2011 Time-averaged Spectra of 4U 1954+319 


Model 

W H 

(10 22 cm -2 ) 

r 

£cut 

(keV) 

±fold 

(keV) 

kT e 
(keV) 

kT 

(keV) 

r 

/ 2 /dof (Prob.) 

2012 A07 ToO observation 

bbody+cutof fpl 
bbody+pow*highecut 
bbody+NPEX 
compTT 

5.4 ±0.1 
4.9 ±0.1 
4.9 ±0.1 
1.7 ±0.2 

1.21 ±0.01 
1.24 ±0.01 
0.77 ± 0.01 

4 23 +0 - 17 
^' z — 0.19 

6.6 ±0.1 

19.9 ±0.5 

21 0 +0 - 9 
Zi,U -0.8 

7 5 +0 - 2 

0.13 ±0.01 
0.13 ±0.01 
0.14 ±0.01 
1.00 ±0.01 

0 q+0.13 
5 -^-0.14 

971.5/696(2.0 x 10" 11 ) 
890.1/696(6.9 x 10" 7 ) 
868.8/695(7.2 x 10“ 6 ) 
780.8/697(0.015) 

2011 A06 observation 

bbodyrad 
diskbb 
cutof fpl 
pow*highecut 
compTT 

1.5 ±0.1 

3.0 ±0.1 

2 1 +0 - 2 
Z '-0.1 

2.0 ± 0.2 

1.6 ±0.1 

-0.73 ±0.17 

rj 90 +O.I 9 

U.Z5 — 020 

3.0 ±0.1 

1.9 ±0.1 
3.0 ±0.3 

<5.1 

1.28 ±0.01 
2.19 ±0.03 

oc+0.10 
U - 8:5 -0.04 

15.3^8 

162.7/124(0.01) 

162.3/124(0.12) 

133.9/123(0.23) 

132.4/122(0.24) 

129.0/122(0.31) 


Notes. All the spectral uncertainties are given at lcr confidence level for single parameters. The photo-absorption and Gaussian components are included in all 
the spectral models. The parameter E cut in the xspec cutoff pi model is actually a folding energy, and shown in the £f 0 id column to compare the highecut 
model. 
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Energy (keV) 

Figure 7. Pulsed fraction of 4U 1954+319 in 201 1 (blue) and 2012 (red). 
(A color version of this figure is available in the online journal.) 



1 10 
Energy (keV) 


Figure 8. (a) Time- averaged, background- subtracted raw XIS-FI, XIS-BI, and 
HXD-PIN spectra of 4U 1954+319 during the two Suzaku pointings in count 
rate space. The modeled NXB, CXB, and 1% of NXB (a typical uncertainty of 
the NXB modeling) are also indicated, (b) Spectral ratio of the 2012 ToO data 
to those from the 2011 observation. The gray diamond of the HXD-PIN point 
represents an error region including 0.3% systematic uncertainties. 

(A color version of this figure is available in the online journal.) 

of 2.1 x 10 -4 ). Then, when we try various abundance tables 
(e.g., angr; Grevesse & Anders 1989; aspl, Asplund et al. 
2009) with the phabs absorption model (Balucinska-Church & 
McCammon 1992), the same comptt continuum changes 
Ah up to by ~50%. Under these uncertainties, we assume 



1 5 10 50 

Energy (keV) 


Figure 9. Ratios of the best-fit spectral model to the 2012 Suzaku data, using 
(a) a single power law, (b) a cutoff power law, (c) a high energy cutoff, 
(d) the NPEX, and (e) the comptt. The 6.4 keV Gaussian emission is included 
in panels (b)-(e), while a soft blackbody component ( kT ~ 0.13 keV) is added 
in panels (b)-(d). 

(A color version of this figure is available in the online journal.) 



Energy (keV) 

Figure 10. (a) X-ray spectra of 4U 1954+319 in vF v form observed dur- 
ing the 2011 quiescent and 2012 flare-up, utilizing the comptt continuum 
(Model C). For comparison, INTEGRAL spectra, JEM-X (3-30 keV) and 
ISGRI (20-80 keV) for the flaring episode in 2008, are overlaid with a comptt 
model fit as used by Marcu et al. (2011), (b) residuals of the comptt model in 
2012. (c) The same as panel (b) but for the 2011 data. 

(A color version of this figure is available in the online journal.) 
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]Vh = 1.7x 10 22 cm -2 , derived from the best- fit comptt model, 
as the default value in the discussion. 14 Subtracting the 21 cm 
interstellar column density in the direction to 4U 1954+319, 
A H = 0.89 x 10 22 cm -2 (Kalberla et al. 2005), the intrinsic 
column density in the binary system is Ah ~ 0.81 x 10 22 cm -2 . 

We also search the HXD-PIN data for a possible presence of a 
cyclotron resonance scattering feature (CRSF), which provides 
a direct probe of the 5-field at the accretion column (Makishima 
et al. 1999; Schonherr et al. 2007). Although a shallow structure 
at —40 keV is visible, the depth of the absorption (cyclabs) 
was constrained to D < 1.1 (lcr). Thus, no CRSF was found in 
the Suzaku data. 

3.5. X-Ray Spectrum in the 2011 Quiescent State 

In the first Suzaku observation performed in 2011 before 
the flare, the source was fainter by more than an order of 
magnitude, with the absorbed 2-10 keV X-ray flux 5.36 zb 

0.03 x 10 -12 erg s _1 cm -2 (comptt model). Since the HXD-PIN 
detection is marginal, we only took the XIS-FI and BI spectra 
into account for this data set. Table 2 summarizes the spectral 
models used to fit the data. Simple power-law and blackbody 
models failed with Xv > 1.3, while the diskbb, cutoff pi, 
pow*highecut, and compTT models were all acceptable. Using 
the same comptt model, the unabsorbed 1-10 keV flux was 
6.38 zb 0.09 x 10“ 12 erg s _1 cm -2 . The lcr upper-limit intensity 
of a narrow (cr = 10 eV) Fe-Ka line at 6.4 keV was 
/ Fe < 9.8 x 10 -7 photon cm -2 s -1 . 

3.6. Luminosity and Mass Accretion Rate 

As seen above, between the quiescent phase in 2011 and 
the outburst phase of 2012 the time-averaged X-ray flux of 
4U 1954+319 changed by a factor of about 100. Assuming 
a distance of 1.7 kpc, the fluxes observed in 2011 and 2012 
correspond to 1-70 keV luminosities of 0.022L 3 5 and 2.IL35 
with L35 = 10 35 ergs -1 (d/1.7kpc) 2 . Correcting for the X-ray 
absorption increases these luminosities to 0.025L 3 5 and 2.2L 3 5, 
respectively. These values fall in the luminosity range typically 
seen in SyXBs, L x — 10 33 -10 35 erg s _1 (Lii et al. 2012). 

We can estimate the accretion rate onto the NS as 

Rm , (4.7 x 10 13 g s -1 in 2011 

i)GM ns 1 3.8 x 10 15 gs _1 in 2012 

where 77 = 0.3 is the efficiency of accretion, 5 ns = 10 km is 
the NS radius, and M NS = 1.4 M© is the NS mass. 

Since 4U 1954+319 is a highly variable source (Section 3.2), 
we further investigated, in Figure 11, distributions of count 
rates in the XIS and HXD-PIN light curves, referring to Vela 
X-l studies by Fiirst et al. (2010). The count rate distributions 
closely follow a log-normal distribution in both the quiescent 
phase (201 1, Figure 1 1(a)) and the outburst (2012, Figures 1 1(a) 
and (b)). Assuming the 2012 time-averaged spectral shape, we 
further converted these to luminosity at the top of Figure 1 1(a). 
As previously shown in Figure 8, the spectrum becomes harder 
as the X-ray luminosity increases from the quiescent phase 
(2011) to the outburst (2012). Such a hardening trend is also 
detected within the 2012 observation when we divide the data 
into three intensity- sorted spectra (Figure 12). 

14 If the Fe-Ka line is due to fluorescence from an optically thin gas which 
spherically surrounds the source, its equivalent width is expected to be linearly 
correlated with the column density, Ah as EW = 100 (Ah/ 10 23 cm -2 ) (Inoue 
1985). The line could therefore be produced in a gas with a column density of 
A h = 2.9 x 10 22 cm -2 . 


4. DISCUSSION 
4.1. Summary of the Observations 

Combining the Swift/ BAT and RXTE/ ASM long-term mon- 
itoring with the Suzaku observations, we derived the following 
observational results which need to be explained. 

1. The slow 5.4 hr NS spin period of 4U 1954+319 exhibits a 
—7% fluctuation over —8 yr with at least four reversals in 
the sign of the period derivative. 

2. In the 2011 October and 2012 November Suzaku observa- 
tions, before and after the 2012 November X-ray outburst, 
the 1-70 keV X-ray luminosities were 2.2 x 10 33 ergs _1 
and 2.1 x 10 35 erg s _1 , respectively. In the former and the 
latter observations, the source was in a spin-down and a 
spin-up phase, respectively. 

3. The 0.5-70 keV X-ray pulse profile shows a main peak 
and a sub peak. The PF mostly increases with energy from 
-60% at 3 keV to 80% at 30 keV. 

4. During the 2012 outburst, the X-ray flux was highly variable 
on 10-1000 s timescales, with many irregular short flares 
each with a typical duration of <100 s. The Fourier 
power density spectrum is dominated by red noise, i.e., 
became redder from —0.77 in 2011 to —1.66 in 2012. 
The X-ray count rates follow a log-normal distribution in 
both observations. 

5. The 0.8-70 keV phase-averaged spectra in the 2012 out- 
burst can be represented by the Comptonization spec- 
tral models; compTT, BB+NPEX, or BB+highecut. A flu- 
orescent Fe-Ka line at 6.37 keV is detected with an 
EW of EW = 28.6 =b 0.3 eV, and a line flux of / Fe = 
7.5 + q 8 7 x 10 -5 photons cm -2 s -1 . The line is narrow, with 
a lcr upper limit for the Gaussian width of 50 eV. The 
quiescent 2011 spectrum is also well described by comptt 
with kT — 0.85 keV, or also by diskbb, cutoff pi, and 
highecut models. 

6. The 2012 XIS spectrum is harder than that measured in 
2011. The intensity- sorted X-ray spectra in 2012 slightly 
hardens with source brightness. 

7. The 1-70 keV band does not contain any statistically 
significant cyclotron lines. 

4.2. Are SyXBs a Subclass ofLMXB? 

X-ray binaries hosting NSs are conventionally classified 
based on stellar types of the optical counterparts. Commonly, 
HMXBs host a strongly magnetized NS and a high-mass com- 
panion (i.e., Be stars or OB supergiants), while LMXBs are 
thought to host a weakly magnetized NS and a low-mass star 
(e.g., K-type stars). However, this old optical identification ig- 
nores accretor properties (e.g., NS 5-field) which mainly de- 
termine the X-ray radiation. In fact, the conventional stereotype 
assuming weakly magnetized NSs in LMXBs is now challenged 
by recent observational evidence of strongly magnetized NSs in 
LMXBs, e.g., 4U 1626-67 (Camero-Arranz et al. 2012) and 
4U 1822-37 (Sasano et al. 2013). 

The required overhaul of this outworn notion is more clearly 
shown in SyXBs. In contrast to the conventional classification 
which makes SyXBs the same category as LMXBs, the apparent 
high PF (60%-80%) of SyXBs (Figure 7, Section 4.1, point 3) 
suggests a similarity to strongly magnetized and highly pulsed 
NSs in HMXBs (PF — 20%-100%; Lutovinov & Tsygankov 
2008) rather than to non-pulsed NSs in LMXBs or to weakly 
pulsed millisecond X-ray pulsars (typical PF at a few percent; 
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Figure 11. Histogram of count rates of the (a) 40 s binned 0.5-10 keV XIS and (b) 80 s binned 15-70 HXD-PIN light curve. Panel (a) shows the summed XIS rates 
in 2011 (blue) and 2012 (red). The corresponding luminosity is shown in the upper axis where the XIS 0.5-10 keV rate is converted to the average 0.5-100 keV 
X-ray luminosity using the 2012 observation and assuming d = 1.7 kpc. The solid curves in both panels show the best-fit Gaussians. Best fit Gaussian mean p and 
standard deviation sigma cr in log-space are derived to be p = —0.31 ± 0.01, a = 0.23 ± 0.01 (2011), p = 1.06 ± 0.01, cr = 3.31 ± 0.01 (2012) for panel (a), 
and p = —0.41 ± 0.02, cr = 0.44 ± 0.01 for panel (b), respectively. The three regimes, bright (A; > 30 counts -1 ), moderate (B; 14-30 counts -1 ), dim phases 
(C; < 14 counts -1 ) and the 2011 quiescent state (D; < 2 counts -1 ) are shaded. 

(A color version of this figure is available in the online journal.) 
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Figure 12. Spectral ratios of 4U 1954+319 to the average one in 2012. Four 
different X-ray luminosity states, (A), (B), (C), and (D) which were defined in 
Figure 11. 

(A color version of this figure is available in the online journal.) 

Patruno & Watts 2012). In addition, if we look at the location 
of SyXBs in the “Corbet diagram” (Corbet 1986) shown in 
Figure 13, SyXBs again resemble HMXBs due to their wide- 
orbit and slow NS rotation. Finally, SyXB X-ray spectra and 
that of 4U 1954+319, in particular, and the other behavior 
summarized in Section 4.1 also suggest a closer similarity of 
4U 1954+319 to X-ray pulsars in HMXBs with a 10 12 -10 13 G 
5-field. 


Even though many of these points argue toward a closer 
similarity of SyXB and HMXB, such an association would 
be a significant challenge to our current understanding of 
X-ray binary evolution. Stellar evolution theory shows that 
the evolutionary timescale for the companion of the NS to 
become an M4-5 III counterpart is much longer than the typical 
lifetime of binary systems with OB-type donors, and thus, the 4U 
1954+319 system would be old. If one dogmatically expected 
NS 5-field decay, the high 5-field indicated by the strong X-ray 
pulsation easily causes a contradiction to the assumption of an 
old binary system. Furthermore, it is also difficult to keep such a 
binary system after a supernovae explosion, since the common 
envelope hypothesis, which is usually assumed to make close 
binary systems hosting NSs in LMXBs, cannot be reconciled 
with the suggested large orbit of SyXBs (Figure 13). Te 
conventional evolutional theory leaves us with a real conundrum. 

In addition to needing to study the population synthesis of 
such binaries (see, e.g., Lli et al. 2012; Chakrabarty & Roche 
1997), the observed nature of SyXB requires us to investigate 
more exotic scenarios, e.g., a magnetar-descendent was captured 
by an evolved M-type giant in a close encounter via a magnetic 
braking interaction between the NS strong 5-field and a large 
atmosphere of the M-type giant. An alternative possibility is an 
accretion-induced collapse of a white dwarf (Nomoto & Kondo 
1991) which potentially produces not only weak- field pulsars 
but also a magnetar (Thompson & Duncan 1995). 

In the following, we take the features of 4U 1954+319 listed 
in Section 4.1 as observational clues to the mystery of SyXB, 
comparing the conventional accretion model and the revisited 
quasi- spherical accretion model, regardless of the conventional 
HMXB and LMXB classification and population synthesis. 

4.3. The X-Ray Spectrum of4U 1954+319 

The first clues on the nature of 4U 1954+319 come 
from broadband spectroscopy. We assume that the observed 
X-ray spectrum is uniquely characterized mainly by M, 
5-field, and inclination of the compact object regardless of the 
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Figure 13. Known SyXBs (red stars) on the Corbet diagram. Sources with either an orbital or rotation period detected are indicated as (red) arrows. HMXBs hosting a 
Be star (green filled circles), OB supergiant (blue filled triangles), and LMXBs (orange open squares; only X-ray pulsars) are also shown from the catalogs (Liu et al. 
2006, 2007; Bodaghee et al. 2007). Arrows on the bottom axis represent atoll and Z-type LMXB systems for which only orbital periods were measured (Liu et al. 
2007). 

(A color version of this figure is available in the online journal.) 


type of the mass-donor. As shown in Section 2.3 and previ- 
ously suggested by Masetti et al. (2007b), spectral modeling im- 
plies that Comptonization is the dominant process forming the 
X-ray continuum, similar to the accretion column of the canon- 
ical X-ray pulsars in HMXBs. 15 Thus, in Figure 14, we now 
compare the X-ray spectrum of 4U 1954+319 with different 
X-ray spectral shapes of other related X-ray pulsars. 

Figure 14(a) shows a comparison of the X-ray spectrum of 
4U 1954+319 with all archival Suzaku SyXB observations: GX 
1+4, 4U 1700+24 (Nagae et al. 2008), and IGR J16194-2810 
(Kitamura et al. 2014). The rotational, orbital, and optical 
information of these X-ray sources is listed in Table 3 with 
the Suzaku observation records and references. To extract these 
spectra, we use the same procedure as for 4U 1954+319 
(Section 2.2). The vF v spectra were obtained by fitting the data 
with a cut off pi model with interstellar absorption (phabs) 
and, if needed, an iron line. GX 1+4 exhibits harder spectra 
peaking at ~20-40 keV than 4U 1954+319. The spectral shape 
of IGR J16194— 2810, for which orbital and pulse periods have 
not yet been detected, is similar to that of 4U 1954+319 in 
quiescence, and thus, their accretion environment is expected to 
be similar. 

In many accreting NSs with 5-fields in the 10 12 G regime, cy- 
clotron lines at energies E cyc = 11.6(5 cyc /10 12 G)(l +z) -1 keV 
have been seen, where z ~ 0.2 is the gravitational red shift 
of the NS (Caballero & Wilms 2012; Schonherr et al. 2007). 
In Figure 14(b), we compare 4U 1954+319 with well-studied 
CRSF sources. The 1.24 s X-ray pulsar Her X-l is a prototyp- 
ical CRSF source with E cyc ~ 36 keV (5 ~ 3.1 x 10 12 G) 
with an unusual optical companion (A9-B companion with a 
mass of M c ~ 2.3 M 0 ). The HMXBs 1A 1118—616 and A 
0535+262 exhibit the strongest 5-fields among known CRSF 


15 Note that previous observations also included a soft excess which could be 
described with plasma emission with kT ~ 50 eV (Masetti et al. 2007b), 
which is not seen here since the Suzaku low energy data below ~1 keV have 
uncertainties due to the contamination on the optical blocking filter. 


sources, reaching B ~ 4.7 x 10 12 G and ~4.1 x 10 12 G, respec- 
tively. IGR J 16393— 4643 was originally classified as a SyXB 
but was finally identified as a HMXB (Bodaghee et al. 2012; 
Pearlman et al. 2011). As shown in Figure 14(b), the spec- 
tral shape of 4U 1954+319 is similar to 1A 1118—616 and A 
0535+262 rather than Her X-l. Although a widely accepted 
physical model has not yet been established to explain the 
X-ray continuum, it has been suggested that due to cyclotron 
cooling the high energy cutoff appears around the CRSF energy. 
This means that the X-ray continuum is expected to extend to 
higher energies as the 5-field becomes stronger (Makishima 
et al. 1999). 

The non-detection of a CRSF in 4U 1954+319 with the XIS 
and HXD means that either the 5-field is not in the range 
5 = 9.0 x 10 10 -6.0 x 10 12 G, or that the cyclotron line is 
too faint to be detected in the data. The latter is not too unlikely, 
with a number of sources with similar behavior being observed. 
A possible explanation for this is that CRSFs can be smeared out 
and/or filled by photon spawning, i.e., the emission of photons 
close to the resonance energy during their de-excitation from 
higher Landau levels (Schonherr et al. 2007; Fiirst et al. 2011a). 
If, on the other hand, 5 > 10 13 G, according to Makishima et al. 
(1999) we would expect a harder X-ray continuum than that of 
HMXBs with a cyclotron line. 

In Figure 14(c), we compare the spectrum of 4U 1954+319 
with those of some HMXBs, which have not shown a clear 
CRSF: HMXB 4U 1909+07, 4U 2206+54, and 4U 0114+65. 
The phase-averaged spectral shapes of these sources are quite 
similar to that of 4U 1954+319. From a spectral view point, they 
can be regarded as triplets. 4U 1909+07 is a HMXB with an OB 
companion hosting a 604 s X-ray pulsar in a 4.4 day orbit. Its 
pulse period is reported to show a random walk-like behavior 
(Fiirst et al. 2011a), which also resembles the period history of 
4U 1954+319 and possibly suggests wind type accretion rather 
than the stable accretion disk. The other two X-ray pulsars, 4U 
2206+54 and 4U 0114+65, exhibit very long rotational periods 
of ~5500 s and ~9700 s, respectively (Finley et al. 1992; 
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Figure 14. Spectral comparison of 4U 1954+319 with different categories of NS binaries: (a) other SyXBs, (b) CRSF sources, (c) long period pulsars, and (d) the 
LMXB Aquila X-l (Sakurai et al. 2012) in the low/hard state. 

(A color version of this figure is available in the online journal.) 


Table 3 

Archival Suzaku X-Ray Pulsars Comapred with 4U 1954+319 


X-ray source 

A 

P orb 

Optical 

ObsID 

Exp 

Note 


(sec) 

(day) 

Type 


(ks) 


4U 1954+31 

18300 

?? 

M4 III 

(Table 1) 

60.2/60.8 

this work 

GX 1+4 

150 

1161 

M5 III 

405077010 

97.2 

cutoffpl [1] 

4U 1700+24 

?? 

(404?) [2] 

M2 III 

402023010 

50.3 

cutoffpl [3] 

IGR J16194-2810 

?? 

?? 

M2 III 

403024010 

45.6 

cutoffpl [4] 

Her X-l 

1.24 

1.70 

A9-B 

100035010 

36.9 

cutoffpl+cutoffpl [5] 

1A 1118-616 

408 

24 

09.5 

403050010 

22.0 

cutoffpl 

A 0535+262 

103 

111 

09.7 

404055010 

31.8 

highecut 

IGR J 16393-4643 

912 

50.2 

B? [6] [7] 

404056010 

50.6 

cutoffpl [8] 

4U 01 14+65 

9720 

11.63 

Blla 

406017010 

106.6 

cutoffpl+cutoffpl 

4U 2206+54 

5555 

19.25 

09.5V 

402069010 

51.6 

cutoffpl 

4U 1909+07 

604 

4.4 

OB 

405073010 

29.2 

cutoffpl 

Aql X-l 


19.0 

K1 

402053030 

19.7 

bbody+compPS [9] 


Notes. The photo-absorption is included in all the spectral models. Applied spectral models and references are noted. [1] Chakrabarty & Roche 1997; 
[2] 404 day period in 4U 1700+24 was not confirmed in Corbet et al. 2008; [3] Masetti et al. 2006b; Nagae et al. 2008; [4] Masetti et al. 2007a; Kitamura 
et al. 2014; [5] Enoto et al. 2008; [6] Bodaghee et al. 2012; [7] Pearlman et al. 2011; [8] Nespoli et al. 2010; [9] Sakurai et al. 2012. 


Reig et al. 2009; Masetti et al. 2006a; Farrell et al. 2008). Thus, 
we expect that such a long rotational period is strongly related 
to the spectral shape in Figure 14(c) of these sources. 

Let us finally compare, in Figure 14(d), the spectrum with the 
low-luminosity (Lx < 10 36 ergs _1 ) low/hard state of LMXBs 
as motivated by Kitamura et al. (2014) who tried to search for 
similarities between SyXBs and low- field LMXB -NS s (e.g., 
Aql X-l) rather than high-field NSs, based on the optical 


classification. Although the L-field of such LMXBs starts to 
dominate the accretion flow and an optically thin spherical flow 
would realize a similar geometry as the SyXB, as we already 
discussed, SyXBs should be different from other LMXBs due to 
(1) slow pulsation with a large PF (Figure 6) suggesting a high 
L-field; (2) the large orbital period (Figure 13) suggesting wind- 
type accretion; (3) the lack of clear evidence on the accretion 
disk; (4) no observations of LMXB -like strong soft blackbody 
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spectra; and (5) the strong narrow 6.4 keV iron fluorescence 
line, which is more typically found in HMXBs but rare (~10%) 
in standard LMXBs (Torrejon et al. 2010). In addition, LMXBs 
are usually characterized by highly ionized 6. 6-6. 9 keV lines. 

In conclusion, the broadband 4U 1954+319 spectrum best 
resembles that of X-ray pulsars in HMXBs, especially slowly 
rotating pulsars (e.g., 4U 0114+65 or 4U1909+07) without an 
apparent accretion disk as well as another quiescent SyXB, 
IGR J16194— 2810. Therefore, a strong field >10 12 or higher 
is favored rather than a low field as in non-pulsating or weakly 
pulsed NSs in LMXBs. 

4.4. The Long Spin Period of4U 1954+319 

Let us now investigate the timing behavior of the source. As 
already discussed previously (e.g., Marcu et al. 2011; Corbet 
et al. 2008, and references therein), the 5.4 hr period of 4U 
1954+319 can only be explained as being the rotation period of 
the NS. Both an orbital modulation or an M-type star pulsation, 
can be ruled out. The long-term pulse history implies at least 
four spin-torque transitions between a spin-up and a spin-down 
phase. The period evolution indicates that outside the outbursts 
the NS is gradually slowing down, while during the outbursts it 
enters a spin-up trend correlated with high accretion rates. This 
behavior suggests that it is reasonable to assume that the source 
is near an equilibrium period of accretion momentum transfer. 

In the standard disk accretion model (Ghosh & Lamb 1979), 
the equilibrium period is Equation (76) in Shakura et al. (2012) 

P eq » s, (5) 

where = M/(10 16 gs -1 ) is the mass accretion rate, and 
/r 3 o = /z/(10 30 Gem 3 ) is the NSs magnetic dipole moment. 
Inserting the values for 4U 1954+319 measured in 2012, 
i.e., Mi6 = 0.38 and 5 e q ~ 19400 s (~5.4 hr) yields an 
extremely high, magnetar-like 5-field reaching B ~ 10 16 G. 
This evaluation would hold even if considering the long-term 
history, since the average mass accretion rate is close to the 
value in the 2012 Suzaku observation. This is clearly shown in 
Figure 15, where the long-term distribution of mass accretion 
rates is illustrated estimated from Figure 2. 

While a magnetar-like 5-field has been suggested for the 
SyXB IGR J16358— 4726 (Patel et al. 2007), the field derived 
above is an order of magnitude higher than the strongest 
magnetar field (SGR 1806—20, with 2.1 x 10 15 G; Kouveliotou 
et al. 1998; Enoto et al. 2010b). 4U 1954+319 also does not 
show any other magnetar-like properties, such as short bursts 
or giant flares (Enoto et al. 2009). In addition, it is unclear 
whether standard disk theory can be applied to SyXBs. For 
example, Gonzalez-Galan et al. (2012) show that the spin 
and X-ray intensity behavior of the SyXB GX 1+4 cannot be 
explained either by standard disk accretion theory or even by 
more advanced ideas such as accretion from a retrograde disk. 
We therefore have to also consider alternative accretion models. 

As the Corbet diagram (Figure 13) implies large orbital 
periods of SyXBs (e.g., 5 0 rb = 1161 d for GX 1+4; Hinkle 
et al. 2006), it is reasonable to assume that 4U 1954+319 does 
not accrete via Roche lobe overflow but is a wind-fed accretor. 
The large orbital period leads to a slow NS orbital velocity. 16 


16 Although the orbital radius of 4U 1954+319 has not yet been detected, 
assuming a circular orbit with its orbital period at 500 days, we will derive the 
semimajor axis at 3 x 10 13 cm via the Kepler’s la w and m asses of the NS and 
M giant. Then the orbital velocity becomes v = y/GM/a ~ 30 km s -1 . 
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Figure 15. Histogram of the long-term 4U 1954+319 X-ray intensity in mCrab 
unit derived from Figure 2. Corresponding X-ray luminosity and mass accretion 
rate are shown in the upper axes, assuming the distance at 1.7 kpc and rj = 0.3. 
See Table 1 for the conversion factor from mCrab to X-ray flux. The two Suzaku 
observations, one INTEGRAL observation (Marcu et al. 2011), and the average 
of the distributions are shown as arrows. 

(A color version of this figure is available in the online journal.) 

Combining this with a slow stellar wind compared to HMXBs 
and with a low mass accretion rate, small angular momentum 
differences are an intrinsic importance of the SyXBs which 
makes the mass accretion more quasi- spherical. 

As recently shown by Shakura et al. (2013), the X-ray 
luminosity of 4U 1954+319 actually puts the source in the 
regime of quasi- spherical, subsonic setting accretion in the 
magnetosphere. Figure 16 shows a sketch of this alternative wind 
accretion scenario for 4U 1954+319, which is possible for 
Lx ^ L th = 4 x 10 36 ergs 1 and M\^ ^ 4. Similar to the 
standard Bondi-Hoyle-Lyttleton accretion, in this scenario, a 
bow shock is formed in the wind at a characteristic distance (5p) 
from the NS, and matter entering this region is accreted onto 
the compact object. Unlike the classical Bondi-Hoyle-Lyttleton 
accretion at higher luminosities, for the low Lx of 4U 1954+319, 
Compton cooling becomes ineffective, and the matter can 
settle in the magnetosphere with a subsonic velocity through 
a quasi- spherical hot shell (Section 1). If this happens, the 
angular momentum can be transferred via large-scale convective 
motions or turbulence in this quasi-static shell. This results in 
an equilibrium spin period of (Postnov et al. 2012) 

^ 13.000 © 

where = u/(1000 kms -1 ) is the relative velocity between 
the stellar wind and the NS. 

The most uncertain parameter in Equation (6) is v$. The 
wind velocity strongly affects 5 eq , but it is not well known for 
the late-type M giants due to lack of observations. Existing 
measurements for wind speeds are mainly available for massive 
OB or Wolf-Rayet stars with high mass-loss rates = 

10 -1 °-10 -5 Mq yr -1 ). These radiatively driven winds are very 
fast (up to ~2500kms -1 ; Lamers & Cassinelli 1999). The 
situation is different for SyXB systems. Here no strong evidence 
for fast winds exists. Late F-M stars are usually thought to 
have smaller mass-loss rates (<10 _1 ° M 0 yr _1 ) with much 
slower terminal velocities (<100 kms -1 ; Espey & Crowley 
2008). These slower speeds are consistent with the lower escape 
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/ I 


Figure 16. Schematic view of the quasi- spherical accretion and subsonic settling accretion (see details in Shakura et al. 2014). Ra and Rb are the Alfven and Bondi 
radii, respectively. 

(A color version of this figure is available in the online journal.) 


velocities of the stars. For example, if the donor in 4U 1954+31 
has 1.2 M q and a radius of 80 5 0 as suggested by Mattana 
et al. (2006), then the escape velocity from the companion star 
would be only ~76kms _1 . To our knowledge, the only direct 
measurement of the wind speed in a SyXB was performed for 
the M6-giant V2116 Oph, the donor star of the SyXB GX 1+4, 
where the wind speed was found to be ~ 100-200 km s -1 
(Chakrabarty et al. 1998; Hinkle et al. 2006). 

Assuming that the wind in 4U 1954+319 is comparable to that 
in GX 1+4, Figure 17 shows the dependence of the equilibrium 
period for several plausible different combinations of F orb , 
M, and B as calculated with Equation (6). If 4U 1954+319 
has a typical accreting NS with a 5-field of 10 12-13 G, then its 
equilibrium period can only be explained by a fast M-star wind 
(e.g., ~300 km s -1 , probably even faster than that observed in 
GX 1+4) and the long orbital period (e.g., ~ 5 00 days). 


A possible explanation for such shots in 4U 1954+319 are 
individual blobby structures in the accretion flow. To see whether 
this interpretation is reasonable, we need to take a look at the 
typical size and timescales of this system (see Figure 17 right 
panel for a visualization). In the model by Postnov et al. (2012), 
the stellar wind is gravitationally captured by the NS inside the 
Bondi radius, 

R b = 2GMNS = 3.78 x 10 10 i^ 2 cm. (7) 

V 2 

For vs = 0.3, i.e., a stellar wind velocity of 300kms _1 , the 
Bondi radius R B ~ 4.2 x 10 11 cm. Inside this radius a quasi- 
static shell is formed above the magnetosphere. This is the 
reservoir of matter from which the NS accretes. In this model 
the co-rotation radius is (Postnov et al. 2012) 


4.5. Timescales of Irregular Flares 

A further handle on the accretion mechanism at work in 
4U 1954+319 comes from the short-term variability of the 
source. The quiescent state in 2011 shows a steady weak 
accretion (in Figures 4(a) and (b)), while the outburst in 2012 
exhibits many irregular flare on timescales of <100 s (Figures 3 
and 4 (c) and (d)). This highly variable fluctuation indicates a 
blobby mass accretion (Figure 16). The X-ray light curves show 
that the X-ray variation is mainly due to short-timescale spiky 
flares. Similar spike-like structures, also called “shots,” were 
studied previously in analyses of the light curves of accreting 
black holes such as Cygnus X-l (Negoro et al. 1994; Yamada 
et al. 2013). In these objects, the spectrum gradually softens 
during the peak of a shot, followed by an immediate hardening. 
Figure 18 shows a similar analysis for 4U 1954+319. The shots 
in the light curves are represented by a shape with a nearly 
symmetric rise and decay (Figures 18(a)-(c)) on a timescale of 
~50 s (at which the count rate drops below the average) with 
signs of a slight hardening in the early phase of the shot-like 
behavior which is indicated in (Figures 18(d) and (e)). 


R c = 


(gmp 2 y /3 

v ~^ r ) 


= 1.21 x 10 11 cm 


5.40 hr 


2/3 


( 8 ) 


Even though the co-rotation radius R c is very large, it is 
still inside the Bondi radius if the stellar wind is slower than 
~600kms -1 (Figure 17 right). This sets an upper limit for the 
wind speed in the system. 

The size of the magnetosphere, characterized by the Alfven 
radius where the gas pressure becomes equal to the magnetic 
pressure (Ghosh & Lamb 1979), is modified for quasi- spherical 
accretion at low mass accretion rates and becomes (Postnov 
et al. 2012; Shakura et al. 2013) 


R a ~ 1.6 X 10 9 cm(^-) . (9) 

\Mi 6 / 

With Mie = 0.38, we find R A = 1.9 x 10 9 cm for B = 10 12 G, 
and 6.7 x 10 9 cm for B = 10 13 G, respectively. This means that 
the Alfven radius is smaller by 1-2 orders of magnitudes than 
R c and 5 B . 
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Figure 17. Left: the NS equilibrium period 5 cq as a function of the magnetic field B and different stellar wind parameters (velocity, mass accretion rate, and fixed 5 0 rb 
at 500 days) according to the quasi-spherical accretion model (Postnov et al. 2012, see Equation (6)). The green box in the background shows the observed equilibrium 
period of 4U 1954+319, 5.0-5. 8 hr. The magenta line at high 5-fields shows the predicted equilibrium period for a standard disk (Equation (5) and Ghosh & Lamb 
1979). Right: the Alfven radius from Equation (9) as a function of the 5-field strength and Bondi radii (Equation (7)) for different wind speeds (1000, 500, and 
300 km s -1 ), assuming the observed M i6 = 1.6. The green box shows the co-rotation radius. The right hand axis shows the free-fall time onto the neutron star from 
each radius. 


(A color version of this figure is available in the online journal.) 


This estimate for the size of the Alfven radius is roughly 
consistent with that found from spectral analysis. While the 
2012 spectrum is Comptonized, during the 2011 quiescent 
state Comptonization does not play a strong role. Approx- 
imating the quiescent spectrum by a blackbody and using 
the Stefan-Boltzmann law, the hot spot size on the NS is 
R hot = 77 m • (d/ 1.7 kpc). Since the hot spot size is deter- 
mined by the location where matter couples onto the 5-field, 
we can use this measurement to estimate the size of the Alfven 
radius (Lamb et al. 1973), 

R a < 10 5 km ( ~ 5.9 x 10 9 cm. (10) 

(lOOm/ 

The measured hot spot radius thus yields a 5-field in the 
10 12 -10 13 G range (Figure 17). 

Depending on the ability of the plasma to enter the magne- 
tosphere, matter from the magnetosphere funnels onto the NS 
poles along the 5-field lines. The free-fall time is given by 

, " = 51 s (to i »ot) ' <M) 

The free-fall timescale from the Alfven radius for 5 ~ 10 13 G 
is ~30 s, close to the typical timescale of individual shots 
(Figure 18). Thus, it is possible to interpret the shot-like flares 
as produced during the free-fall of the blobs from the shell 
above the magnetosphere. Even though the cooling timescale 
of the matter is complicated when considering different cooling 
processes, if we assume Compton cooling as an effective process 
(Shakura et al. 2012), the cooling time, 

Looi = 10 s • —7— ^ 10 ^ (12) 

M 16 V10 iU cm/ 

becomes ~1200 s at 5 = 10 13 G and M \ 6 = 0.38. This is 
longer than the typical free-fall time (f coo i > %), which satisfies 
the assumption of a low-luminosity quasi-spherical accretion. 

The log-normal distribution of the X-ray light curve count 
rates (Figure 11) is also consistent with the picture of quasi- 
spherical magnetospheric accretion. As suggested by Fiirst 


et al. (2010), such a log-normal distribution implies that the 
underlying physical process should be multiplicative rather than 
additive, since addition of many independent emissions will 
produce a Gaussian distribution due to the central limit theorem. 
Such a process could be induced via a positive feedback 
loop, where blob accretion induces subsequent accretion, e.g., 
via enhanced Compton cooling, potentially creating positive 
feedback until the matter at the boundary of the magnetosphere 
is exhausted. Alternatively, the structures leading to a log- 
normal distribution of the accretion rate could also be due to 
a highly inhomogeneous and clumpy stellar wind from the 
M-giant companion or due to turbulent processes in the spherical 
shell itself. 

4.6. Absorption in the Spherical Shell 

Finally, let us compare the above accretion scheme with the 
absorption properties of the circumstellar material. These come 
from both the absorption modeling of the X-ray continuum and 
the analysis of the narrow Fe-Ka line at 6.4 keV. The latter 
was detected with Suzaku in the 2012 observation and has been 
previously reported with BeppoSAX (EW = 5 1 zb 20 eV) and 
RXTE (48+ 13 0 eV; Masetti et al. 2007b). A narrow iron line 
feature is characteristic to HMXB-NSs rather than LMXB-NSs. 

As a plausibility check, we assume that Ah is directly related 
to the properties of the shell in the low-luminosity wind accretion 
model. Using the Suzaku- derived Ah, the number density (nu) 
inside the spherical shell (5 < 5 B ) becomes (nu) = Ah/5 b = 
1.9 x 10 10 cm -3 (i>8 = 0.3), corresponding to a mass density of 
(p) = (nn)m v ~ 3.2 x 10 -14 gcm -3 , where m p is the proton 
mass. If the NS is fed by blobs of matter which have the density 
of the shell at the Alfen radius, then in order to reach the average 
mass accretion rate (M i b ) = 4/3tt Rl loh p/t S h 0t , where 5 b i ob 
and /shot are the typical blob size and shot timescale, we find 
5 b iob ~ 8.0 x 10 9 cm with (M\e) = 0.38 and f S hot ~ 50 s. It 
is intriguing that this initial blob radius is comparable to the 
Alfven radius (e.g., R A = 6.8 x 10 9 cm at 5 = 10 13 G) and is 
much smaller than the quasi-spherical shell size 5 B . Thus, the 
shell can store enough mass to generate the observed multiple 
flares and to sustain accretion for longer time intervals. 
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Figure 18. Background- subtracted stacked shot profiles of 4U 1954+3 1 in (a) the 
0.5-4.0 keV, (b) 4.0-10.0 keV, and (c) 15.0-70.0 keV bands. Shot peaks were 
determined from local maxima within ±200 s duration and with a threshold 
factor of / = 1.8 (Yamada et al. 2013). The profiles in panels (a)-(c) are 
normalized by the average rate in individual energy bands. Panels (d) and (e) 
show the divisions of the profiles in panels (b) and (c) by that in panel (a). 


Finally, as a cross check, we note that the ionization parameter 
of the shell at the Bondi radius is § = L x /{nn)r 2 = 2.1 x 
10 4 (r/10 10 cm) -2 ~ 11, i.e., the outer part of the shell is at most 
moderately ionized, making it possible to produce the 6.4 keV 
line in the shell. 


5. CONCLUSION 

We analyzed two Suzaku data sets, obtained in 201 1 and 2012, 
of the SyXB 4U 1954+319 consisting of the slowly rotating 
5.4 hr X-ray pulsar and a late type primary star (M4-5 III). Our 
main results are as follows. 

1. We reconfirmed the slow rotation period and the large 
period fluctuations. We also identified several additional 
features of 4U 1954+319, including a hard continuum 
without a clear CRSF, a narrow Fe-K a line, and sporadic 
flares. The continuum was harder when the source was 
brighter. 

2. From spectral and timing properties, 4U 1954+319 is 
suggested to have a strong 5-field of >10 12-13 G. This is 


supported by the high PF of 60%-80%, the slow rotation, 
the hard featureless continuum, and the narrow Fe-Ko' line. 

3 . The NS is expected to be close to equilibrium because of the 
four spin transitions and large period fluctuations. Such an 
evolution can be explained by a quasi- spherical accretion 
in a subsonic accretion regime which is considered to 
apply to 4U 1954+319 because of its low luminosity 
(0.023-2. 1 x 10 35 erg s -1 ), presumably wide orbital period, 
and rather low wind velocity from the M-type giant. 

4. Recurrent irregular flares during the outburst have a typical 
timescale of ~50 s, which is interpreted as intermittent 
accretion from the Alfven radius of the ~10 13 G NS. 
The log-normal distribution of the light curve indicates an 
underlying multiplicative process (e.g., positive feedback) 
in the accretion. 

5 . The extreme magnetar-like extreme 5-field ( 1 0 16 G) derived 
from the standard disk accretion in Equation (5) is not 
required if we assume a low luminosity and quasi- spherical 
accretion. 

6. The presence of a system like 4U 1954+319, which is very 
likely to involve a high-field NS, provides another example 
that urges us to challenge the conventional classification of 
NS binaries into LMXBs and HMXBs. 
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